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1. Initialize the loop 

lpl = { 

{LCLi, LCLj) , (LCLi, UCLj) , 
{UCL if UCLj), (UCLi, LCLj)', {LCL if LCLj) 
}, 

k = 1. 

2. Create the half space hs defined by the specification k: 

hs = {X e R n I f k (X) > LSL k ] . 

For two-side specifications, repeat Step 3~5 for C7SL*. 

3. Set lpO = lpl, 

v a = vi, and 
v b = v 2 . 

Empty lpl; .-' 

4. If both v a and v b are inside hs then 

add into lpl; 
else 

if one of v a and v b is inside hs then - ^ 

calculate the intersection point v c of v a v b and hs 
add v c into the new loop lpl 

if v b is inside hs then 
add Vi, into lpl . 

5. Set v a = v a+I and 

Vb = Vb + l . 

If not all vertices are done, go to step 4. 

6. Jt = £ + 1. If k < n, go to step 2. 

7. Add the first vertex into the end of lpl to finish the 
loop. The region surrounded by lpl is the feasible space 
of xi vs. Xj with other parameters fixed. 
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1. Choose n distinctive constraints from all in 
specifications and n parameter limits; 

2. LU decompose F into LXJ; 

3. Constitute one distinctive right-side Y from 
specifications and parameter limits. 

4. Solve LUX = Y by forward and back substitution. If the 
solution X* satisfies all the specifications, add it into 
the extreme point list. 

5. If all distinct Ys have been operated upon, go to next 
step. Otherwise go to step 3. 

6. If all distinct constraint combinations have been 
operated upon, go to next step. Otherwise go to step 1. 

7. All extreme points and the facets are now available. The 
feasible decision space and performance space can be 
acquired by projecting the facets into each 2-dimensional 
viewing pane, though other projections and 
representations may be beneficial. 

8. The boundary of the feasible space in each viewing pane 
is a convex hull of the extreme points. Graham-scan 
algorithm can be adopted to retrieve the convex boundary 
in Q(hlg(h)) time where h is the number of all vertices. 
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